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ABSTRACT 

We find the dispersion relation for tightly wound spiral density waves in the surface of rotating, 
self-gravitating disks in the framework of Modified Gravity (MOG). Also, the Toomre-like stability 
criterion for differentially rotating disks has been derived for both fluid and stellar disks. More specif¬ 
ically, the stability criterion can be expressed in terms of a matter density threshold over which the 
instability occurs. In other words the local stability criterion can be written as Sg < Scrit('Cs) k, a, /ig), 
where Ecrit is a function of Vs (sound speed), k (epicycle frequency) and a and /ig are the free parame¬ 
ters of the theory. In the case of a stellar disk the radial velocity dispersion ar appears in Ecrit instead 
of Vs- We find the exact form of the function Ecrit for both stellar and fluid self-gravitating disks. 
Also, we use a sub-sample of THINGS catalog of spiral galaxies in order to compare the local stability 
criteria. In this perspective, we have compared MOG with Newtonian gravity and investigated the 
possible and detectable differences between these theories. 
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1. INTRODUCTION 

Disk galaxies encompass non-axisymmetric features 
such as spiral arms and bars which have significant in¬ 
fluences on the galactic evolution (Sellwood 2014, Ko- 
rmendy & Kennicutt 2004). These spiral arms trigger 
star formation in the outer part of the disk galaxy, drive 
secular changes in the stellar orbits and change the mass 
distribution within the disk. Understanding the origin 
and evolution of these spiral structures is one of the hard¬ 
est problems in astrophysics. However, it is widely be¬ 
lieved that spiral patterns are gravitationally driven den¬ 
sity waves (e.g., those produced by a spontaneous distur¬ 
bance or, in some cases, a companion system) and prop¬ 
agating in the surface of the galaxies. Therefor the den¬ 
sity wave theory is a indispensable tool in this subject. 
Assuming that density waves are small perturbations on 
the density of the galactic disk, one can use linear anal¬ 
ysis in order to find the linear oscillatory modes of the 
disk. However, analyzing the behavior of a density wave 
in a galactic disk and calculating the relevant oscillatory 
modes is a difficult task. The main reason for this diffi¬ 
culty is that gravity plays the dominant role in the evo¬ 
lution of the density waves. On the other hand, gravity 
is a long range force. Correspondingly, perturbations in 
all parts of the system are coupled. This coupling makes 
the analysis very difficult and consequently the oscilla¬ 
tory modes can be found analytically only for Kalnajs 
disks (Kalnajs 1972). Also numerical mode calculations 
have been performed for some simple models including 
the isochrone disk, the Kuzmin disk, and power-law disks 
(Zang 1976; Kalnajs 1978; Vauterin & Dejonghe 1996; 
Pichon & Cannon 1997; Evans & Read 1998; Jalali & 
Hunter 2005). 

In the early 1960s it was understood that for tightly 
wound spiral density waves (waves whose wavelength is 
much less than the radius), the long-range coupling of 
gravity is negligible, and the response of the background 
disk to this tightly wound spiral density waves can be 
determined locally (Lin & Shu 1964). In fact, in this 
approximation, known as WKB approximation, the rel¬ 


evant analysis is completely analytic and independent of 
the given model for the background disk. In other words, 
using this approximation, one can easily determine the 
stability of the given disk against axisymmetric pertur¬ 
bations (density wave). At an arbitrary location in the 
stellar disk, stability requires that 


Qs 


KOr 

3.36GEo 


> 1 


( 1 ) 


where k is the epicyclic frequency, Ur is the radial veloc¬ 
ity dispersion, Eg is the disk unperturbed density, and 
G is the gravitational constant. The requirement m is 
known as the Toomre’s local stability criterion and has a 
simple physical interpretation. In fact this criterion rep¬ 
resents a competition between the stabilizing influences 
of pressure (tJr) and the angular momentum (k) against 
the destabilizing effect of gravity (Eg). In the case of 
a fluid (gaseous) disk, this criterion takes the following 
form 


Qg — 
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> 1 


( 2 ) 


where Vs is the sound speed. It is worth mentioning 
that it was first shown by Safronov (1960) that a thin 
fluid rotating disc can be unstable against local axisym¬ 
metric disturbances under the effect of its own grav¬ 
ity. These criteria, i.e. m and have been modi¬ 
fied for taking into account various effects. For exam¬ 
ple the effect of the thickness of the disk has been in¬ 
vestigated by Toomre (1964) and Vandervoort (1970). 
Two component (fluid -I- stars) disks of finite thickness 
has been considered by Kato (1972), Berlin & Romeo 
(1988), Romeo (1992), Wang & Silk (1994) and Romeo 
& Wiegert (2011). See also Shadmehri et al. (2012) for 
a two component study where the gaseous component 
is turbulent. Rafikov (2001) studied the general case of 
multiple components. The stability of a two-fluid disk 
was investigated by Jog & Solmon (1984), Elmegreen 
(1995) and Jog (1996). Effects of magnetic field on lin¬ 
ear gravitational instabilities in thin galactic disks have 
been studied by Elmegreen (1987, 1994), Gammie (1996), 




2 


Fan & Lou (1997), Kim and Ostriker (2001). Gammie 
(1996) has also considered the effect of viscosity to the 
local stability. Also the local stability in two-component 
galactic disk with gas dissipation has been studied by 
Elmegreen (2011). 

The purpose of this paper is to introduce a modi¬ 
fied stability criterion in the context of Modified gravity 
(MOG) for a single component galactic disk (fluid or stel¬ 
lar) . This theory is a covariant modification of Einstein’s 
general relativity (GR) and has been introduced for ad¬ 
dressing the dark matter problem (see Moffat 2006). In 
the weak field limit, the Poisson equation in this theory 
is different from that of Newtonian gravity. Therefore, 
one might naturally expect that the local gravitational 
stability criterion is different from ((T|) and ([2]). Eind- 
ing the generalized version of the Toomre’s criterion in 
modified theories of gravity is not a new idea. For ex¬ 
ample, Milgrom (1989) has found the local stability cri¬ 
terion in modified Newtonian dynamics (MOND). Also 
Toomre-like criterion in the context of metric f{R) grav¬ 
ity has been calculated for both gaseous and stellar disks 
by Roshan & Abbassi (2014b). Furthermore, this issue 
has been investigated in the context of some modified 
gravity theories which introduce a Yukawa like term in 
the gravitational force by Habibi et al (2014). In this pa¬ 
per (Habibi et al 2014) the local stability of disks in the 
context of MOG, has been also studied briefly as an spe¬ 
cial example for the above mentioned theories. Also some 
approximate stability criteria has been derived. Here, we 
give a more careful treatment and find the exact form of 
the stability criteria and use the relevant observational 
data to find out a detectable difference between Newto¬ 
nian gravity and MOG. 

2. MODIFIED GRAVITY (MOG) 

In this section we briefly review the current status of 
MOG among other extended theories of gravity. Also, 
we briefly review its weak field limit which is necessary 
to achieve the main goal of this paper. 

As we mentioned in the introduction, MOG is a fully 
relativistic and covariant generalization of GR. Moffat 
(2006) introduced this theory as an alternative to dark 
matter theories. MOG is much more complicated than 
GR in the sense that its associated gravitational fields are 
more than GR. In fact, MOG is a Scalar-tensor-vector 
theory of gravity, while GR is a tensor theory. In other 
words, in GR the components of the metric tensor are 
the gravitational fields of the theory. On the other hand, 
in MOG in addition to the metric tensor, there are three 
scalar fields {uj{x^), and G{x^)) and also a massive 
Proca vector field ( in the current literature of MOG, 
the scalar field w is assumed to be constant). Therefor, 
MOG has equipped with some degrees of freedom in or¬ 
der to handle the dark matter problem without invoking 
dark matter particles. This theory has been successfully 
applied to explain the rotation curves of spiral galaxies 
and the mass discrepancy in the galaxy clusters (Brown- 
stein & Moffat 2006, 2007; Moffat 2006; Brownstein 2009; 
Moffat & Toth 2008, 2009, 2013; Moffat & Rahvar 2013, 
Moffat & Rahvar 2014). 

Eor the details of deriving the weak field limit of MOG, 
we refer the reader to Moffat & Rahvar (2013), Roshan 
& Abbassi (2014a). In the weak field limit of this theory. 


the particle’s equation of motion can be written as 


(Pr 

dP 


= -V$ 


( 3 ) 


where <i> is a effective gravitational potential and is given 
by 

$ = 4' -f (4) 


where x is a coupling constant and 0° is the zeroth com¬ 
ponent (or time component) of the vector field. Further¬ 
more, 4/ and (j)^ satisfy the following equations 


V24' = 4:n{l + a)Gp (5) 

(V^ - = -47raGp (6) 


where po and a are the free parameters of the theory 
in the weak filed limit, G and p are the gravitational 
constant and the matter density, respectively. In fact, 
/io is the background value of the scalar field p and a 
is related to the coupling constant x as a = x^/^oG. 
Note that ujq is the background value of the scalar field 
uj. For moe details see Roshan & Abbassi (2014a). The 
observational values of the free parameters a and po are 
known from rotation curve data of spiral galaxies. It 
has been shown by Moffat & Rahvar (2013) that the 
best values for these parameters are a = 8.89 ± 0.34 
and pq = 0.042 ± 0.004fcpc“^. It is worth mentioning 
that, equations (O and dH) can be combined to write the 
generalized Poisson equation as 

/ p-/io|r-r'| 

-y—^p{v')(Px' (7) 


For considering the dynamics of a self-gravitating fluid 
system in the context of MOG, in addition to the gen¬ 
eralized Poisson equation 0, we need the generalized 
version of the continuity and Euler equations. It can 
be shown that the mathematical form of these equations 
are the same as in Newtonian gravity (Roshan & Abbassi 
2014a). The continuity and Euler equations are: 


dt 


+ V • (pv) = 0 


(9v 


-I- (v • V)v 


Vp 

- - - V4> 


( 8 ) 

( 9 ) 


where p is the pressure and v is the velocity of the fluid. 
Equations 0- ® combined with the equation of state of 
the fluid, makes a complete set of equations for describing 
the dynamics of a fluid system in the framework of MOG. 
Equivalently, one can use equations dH) and ® instead 
of dZl). 

On the other hand, in order to investigate the dynamics 
of a stellar system we need the Boltzmann equation and 
the generalized Poisson equation 0. The Boltzmann 
equation in MOG has been investigated in Roshan & Ab¬ 
bassi (2014a), Moffat & Rahvar (2014). Again, one can 
verify that the Boltzmann equation’s mathematical form 
does not change . The main difference is that the effective 
potential 4> appears in the Boltzmann equation instead 
of the Newtonian gravitational potential Therefore, 
the collisionless Boltzmann equation in MOG is 

df df 

^-kv.V/-V4>.^=0 


( 10 ) 
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where / is the phase-space distribution function. In the 
section m we will use (fTOl) for studying the local stability 
of a stellar disk. 


3. DISPERSION RELATION FOR A 
SELF-GRAVITATING FLUID DISK IN MOG 


In this section we study the behavior of a tightly wound 
density wave propagating in the surface of a fluid disk in 
the context of MOG. First, we use the modified Pois¬ 
son equation © in order to calculate the gravitational 
potential of a tightly wound surface density. Then, by 
linearizing the field equations, we find the relevant dis¬ 
persion relation. We assume that the background disk is 
razor-thin and axisymmetric and its plane corresponds 
to the X — y plane. Also, we assume that the disk has a 
barotropic equation of state as p = KYP, where K and 
7 are constant real parameters. The model assumes a 
non-rotating cylindrical system (r, ip, z), such that z axis 
coincide with the rotation axis of the disk and the angle 
ip increases in the direction of rotation. In the cylindrical 
coordinate, by assuming that the fluid is axisymmetric, 
continuity and two components of Euler equation can be 
written as: 


dvr 


dv^ 




dVr V^p dVr '^ip _ ^ 

'"dr r dp r 


dv. 


dr 


ip ^ Xp dvp ^ VpVf 


dp 


dr 
r dp 


( 11 ) 

+ h) (12) 

($ + /i) (13) 


where Vr, Vp and E are the velocity components in the 
radial and azimuthal directions and the surface density, 
respectively. Furthermore, h is the specific enthalpy de¬ 
fined as h = f dp/Y. Now, we linearize equations (llip - 
(fTSl) by assuming that: Y = Eq-I-Ei, Vr = Vro+Vri = Vri, 
Vp = Vpo + Vpi, $ = $0 -|- $1 and h = ho + hi. The sub¬ 
script ” 0 ” refers to the background value of the given 
quantity and ” 1 ” refers to the corresponding perturbed 
quantity. After linearizing, the fluid equations (HHl-dTl 
can be written as (Binney & Tremaine 2008): 


dYi 


(Eorwri) + 

r or Oif 


Eq dvpi 
r dp 


= 0 
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dp 
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(14) 

(15) 

(16) 


where D(r) is the angular rotation rate and the epicyclic 
frequency n is defined as 


K(r) = + 4112 ( 17 ) 

V dr 

Also, equations ([5]) and (| 6 ]) can be linearized as 

= 47 r(l-ha)GEi 5 (z) (18) 

(V2 - = -AnaGYidiz) (19) 


Q^^i(kr+mp uit) ^ Where w is the oscillation frequency 
and A: = ^ is the wavenumber. For the local stability 
analysis we use the Wentzel-Kramers-Brillouin (WKB) 
approximation (or tight winding approximation) which 
requires that ^ ^ 1 (Binney & Tremaine 2008). This 
approximation allows us to neglect terms proportional 
to 1/r compared with the terms proportional to k. Em¬ 
ploying this approximation, equations (Hl-dlll) can be 
combined to give 

(toD — Uj)Ya + kYoVra = 0 (20) 


'^ra — 


(mfl — a;)fc(4>a -|- ha 


Gipa 


A 

2iB 

-? 

oj — mfl 


( 21 ) 

( 22 ) 


where A and B (the Oort’s constant of rotation) are 


A = — (mil — a;)2 (23) 


In order to complete the local stability analysis, let us 
find the gravitational potential of a tightly wound spiral 
perturbation, Ei = j^^gdkr+mp-uit) neighborhood 

of an arbitrary point {ro,po). Locally, the WKB den¬ 
sity perturbation can be considered as a plane wave with 
wavenumber k = /ce^, where is the unit vector in 
the radial direction. Therefor, finding the potential of a 
WKB perturbation reduces to finding the gravitational 
potential of a plane density wave. In order to find the 
gravitational potential of this perturbation, without loss 
of generality, we choose the x axis to be parallel to k(ro). 
Then we use equations (IT^ and (IT^ to find the gravita¬ 
tional potential. To do so, we guess the solution to (fT^ 
as 

= (25) 


where 4)^ and C S'l'e arbitrary constants. The disk is 
razor-thin and so there is no matter outside {z ^ 0 ) the 
disk. Therefore, at z 7 ^ 0 equation m can be written as 
V2(()5 = Po^i- Substituting the potential (051) into this 
equation, one can easily show that C = ±-\/fc^ + /io- On 
the other hand, since matter is located at z = 0 , deriva¬ 
tive of (pi with respect to z is not continuous on the disk. 
In order to fix the parameter we integrate equation 
(fT^ with respect to z in the interval z = —^ to z = -|-^, 
where ^ is a positive constant, and then let ^ 0. The 

result is 


27rGa 


(26) 


a similar procedure for equation (IT9l) gives 


27rG(l + a)^ 
" |fc| 


(27) 


finally, the gravitational potential of a plane wave in the 
context of MOG is given by 


where 6{z) is the Dirac delta function. We 
choose infinitesimal disturbances of the form Qi = 


$ = 
^n. — 


27rG 

~w 


I -I- a — 


s/k^ + pl) 


(28) 
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Equations O-dUl) and (1^51) combined with ha = 
form a closed set of linear equations. We solve 
these equations to obtain the dispersion relation. Using 
equations (EU and (E51) we can find a purely algebraic 
expression between Vra and E^ as: 


^ra 


{mfl — uj)k 
A 




+ a 


|fc|Q: A 
+ ho j 


(29) 

substituting this equation into (I20L and restricting our¬ 
selves to the axisymmetric (m = 0) perturbations, we 
obtain the following dispersion relation 


w 


2 


K^ + k'^vl-2TTGTsQ\k\ 



\k\a \ 
Vk^ + hlj 


(30) 


We derive the local gravitational stability criterion by 
using this dispersion relation. As expected, in the limit 
/rg —>■ 0 or a —>• 0, the dispersion relation for WKB per¬ 
turbations in Newtonian gravity is recovered. Since all 
the quantities in the right-hand side (RHS) of (1^ are 
real quantities, then the disk is stable against local dis¬ 
turbances if > 0 and unstable if < 0. 

In order to find the local stability criterion, we write 
the dispersion relation (1501) in a dimensionless form. To 
do so, we define the dimensionless wavenumber X, and 
the parameter /3 as 



ho 



(31) 


therefore the dispersion relation can be rewritten in the 
dimensionless form 


= 1 -f 


2/3|A| 

Qg 


(1 -b Of) -b 


2^ 

^Vl + A2 


(32) 


where Qg is the dimensionless Toomre’s parameter given 
by equation ([5]). Therefore, the stability criterion > 0 
takes the following form 


Qg > 


2/3 


1 -b /32A2 


(l + a)|A| 


TT+x^j 


(33) 


This is the exact criterion for the stability of the fluid disk 
against a perturbation mode with given wavenumber X. 
Once we have the observational values of a and /ig, and 
also the magnitude of the wavenumber A, then we can 
explicitly determine the required value of the Toomre’s 
parameter for providing the stability at the given point. 
More specifically, we should know the magnitude of the 
ratio Vs/K at the given location. It is natural that if we 
find the maximum value of the RHS of (l33)l for given a 
and /3, then we will get a criterion which is independent 
of the value of the wavenumber. Since this can not be 
done analytically, we have listed the local stability crite¬ 
rion for various values of a and /3 in Table [T] It is worth 
remembering that, by fitting MOG to observed rotation 
curves of galaxies, the free parameters a and /ig have 
been determined to have the values: a = 8.89 ± 0.34 and 
/rg = 0.042 ± 0.004 kpc~^. So in the Table [p we have 
chosen a as 8.89 ± 0.34. On the other hand, in order to 
choose rational values for /3, we have estimated it in the 
solar neighborhood. In the solar neighborhood the epicy¬ 
cle frequency is k ^ 37 km s“^kpc“^, and the radial ve¬ 


Table 1 

Stability criterion for fluid and stellar disks, for 
different values of a and {3. 


a 

h 

FLUID DISK 

STELLAR DISK 

9.23 

0.02 

Qg > 1.0019 

Qs 

> 1.0014 

9.23 

0.04 

Qg > 1.0075 

Qs 

> 1.0074 

9.23 

0.2 

Qg > 1.4271 

Qs 

> 1.3502 

9.23 

0.4 

Qg > 2.5928 

Qs 

> 2.4840 

8.55 

0.02 

Qg > 1.0017 

Qs 

> 1.0013 

8.55 

0.04 

Qg > 1.0069 

Qs 

> 1.0068 

8.55 

0.2 

Qg > 1.3536 

Qs 

> 1.2826 

8.55 

0.4 

Qg > 2.4459 

Qs 

> 2.3452 


locity dispersion is (Jr ^ 38 km s“^. Therefore /3 ~ 0.04. 
It is obvious from Table [T] that in the solar neighborhood 
there is no significant difference between the local stabil¬ 
ity criterion of MOG and Newtonian gravity. Also, it is 
clear that the difference between these theories consider¬ 
ing the local stability will be significant in the locations 
where the ratio of the sound speed to the epicycle fre¬ 
quency is large. In fact, Vs/k is an important length 
scale for comparing MOG with Newtonian gravity in the 
disk galaxies. If this length scale satisfies the inequality 
Vs/k > 2.4 fcpc, the deviation between the local stability 
criterion of MOG and Newtonian gravity will be large. 
In the end of this section, we will discuss the relevant 
observational data for this length scale. 

It is interesting to plot the dispersion relation for the 
growing (unstable) modes and discuss the growing rates. 
Our purpose is to compare the growth rates in MOG and 
Newtonian gravity. To do so, we rewrite the dispersion 
relation (IHHI) as follows 


M 

Qg 


^I -b a - 


Qkl \ 


(1 + 7^) 


(34) 


where the dimensionless wavenumber q is defined as 
q = kvs/K , and s' = iuj/K. Figure [T] shows the fluid 
dispersion relation for various values of Qg and /3. The 
dotted curves are for Newtonian gravity. The maximum 
growth rate in MOG is higher than that of Newtonian 
gravity, and also occurs at smaller wavenumber. As in 
Newtonian gravity (dotted curves), the growth rates in¬ 
crease with decreasing Qg. It is also clear form Figure [U 
that, in MOG the range of instability is stretched to 
slightly larger wavenumber. Furthermore, by increasing 
the parameter (3 the range of instability and also the 
maximum growth rate increase. 

It is worth mentioning that in metric f{R) gravity, by 
increasing the corresponding parameter /3, the growth 
rate decreases (Roshan & Abbassi 2014b). In this theory 
the parameter /3 can be defined as /3 = moVs/n, where 
mg is one of the free parameters. It is surprising that, 
although both MOG and f{R) gravity lead to a stronger 
gravitational force relative to Newtonian gravity, MOG 
increases the growth rate while f{R) gravity decreases 
it. 

It is also instructive to illustrate the neutral stabil¬ 
ity curves, i.e. w = 0, for axisymmetric tightly wound 
spiral density waves in a fluid disk. These curves deter¬ 
mine the boundary of the stable and unstable perturba¬ 
tions. Let us define a critical wavelength Acrit = 27r/A:crit 
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Figure 1. Solutions to the dispersion relation between the 
squared growth rate and dimensionless wavenumber q for a 
fluid disk. Solid curves correspond to Qg = 0.3 and the dashed 
curves are for Qg = 0.5. Higher value for parameter /3 leads to 
greater instability at high wavelength. 



where kcvit — In Newtonian gravity Acrit 

is the largest unstable wavelength for a fluid disk with 
zero sound speed. Using this definition and the disper¬ 
sion relation ([3Q|l. the line separating stable and unstable 
perturbations is the solution of the following equation 


y 

Figure 2. The boundary of stable and unstable axisymmetric 
perturbations in fluid and stellar disks. The chosen value for a is 
9.23 and (3 varies from 0.02 to 0.2. The dotted line corresponds 
to Newtonian gravity and the dashed curves correspond to a fluid 
disk in MOG (different values of /3). Furthermore, the solid lines 
show the neutral curves of a stellar disk. 


Qg{v) 



(35) 


where y = A/Acrit and A = 27:/k. We have shown the 
neutral curves for various values of j3 in Figure [2] Fur¬ 
thermore, in this figure, we have assumed a = 9.23. The 
dotted curve corresponds to Newtonian gravity. Also, 
the dashed curves show the neutral curves of a fluid disk 
and solid lines are the corresponding lines for a stellar 
disk. It is clear that, at small wavelengths, the neutral 
curves of MOG and Newtonian gravity coincide. This co¬ 
incidence is almost true for the interval Q < y < 0.5, i.e. 
A < 0.5Acrit- However, for y > 0.5 deviations between 
MOG and Newtonian gravity appear. More specifically, 
the range of instability is stretched to larger wavelength. 
It is surprising that the largest unstable wavelength in 
MOG is more than 10 times the Newtonian case. In fact 
one can easily show that the largest unstable mode in 
MOG occurs at A = (1-1- a)Acrit- This is a quiet big 
difference between these theories. However, one should 
note that at large wavelengths, there is a serious doubt 
for the validity of the WKB approximation. For example, 
for solar neighborhood Acrit 10 kpc. This wavelength 
is larger than the solar distance to the center of the disk 
galaxy, r ^ 8.5 kpc . Therefore fccritF ^5.3 and for the 
largest unstable wavelength in MOG we have kr ^ 0.5. 
Thus, obviously, the WKB approximation does not sat¬ 
isfied. Therefore, the wideness of the range of instability 
in MOG compared with the Newtonian gravity, can not 
be considered as a significant and detectable difference. 

On the other hand, it is clear in Figure [5] that increas¬ 
ing the parameter /3, leads to a higher Qg for stability of 
all wavelengths. In other words, in MOG, higher Qg pa¬ 
rameter compared with Newtonian gravity is needed to 
overcome the local gravitational instability. This result 


is also consistent with the Table [TJ As we discussed, this 
difference between MOG and Newtonian gravity could 
lead to a detectable deviation. 

The origin of this difference can be simply understood. 
Using the modified Poisson equation 0, the gravita¬ 
tional force between two point masses mi and m 2 that 
are located at ri and r2 respectively, is given by 

F = Fn+ aFN{l-e-i^o^{l + por)) (36) 

where Fjv = —Gmim 2 /r^er is the Newtonian gravita¬ 
tional force and r = ri — r2 and is the unit vector in 
the direction of r. Taking into account the current val¬ 
ues of a and po, one can easily verify that F is stronger 
than Newton’s gravitational force F^r. Now, keeping in 
mind the physical meaning of the gravitational Jeans in¬ 
stability, it is natural to expect more gas pressure for 
preventing the gravitational collapse. In other words, it 
seems like that, when the gravitational force is stronger 
than the Newtonian gravity, then we need higher Qg pa¬ 
rameter. However, this interpretation is not always true. 
For example, in the case of metric f{R) gravity, where 
the corresponding gravitational force in the weak field 
limit is stronger than Newtonian gravity, the required 
value for Qg is smaller than the Newtonian one. 

4. DISPERSION RELATION FOR A 
SELF-GRAVITATING STELLAR DISK IN MOG 

Stability of stellar disk was first studied by Toomre 
(1964), Kalnajs (1965), and Lin & Shu (1966) who 
pointed out that in certain regimes the stellar disk is 
stable against disturbances. When dealing with stellar 
disk stability, we need to use the collisionless Boltzmann 
equation m, as well as the modified Poisson equation 
0. Since we assume that the stellar disk is collisionless, 
the stars can not create any pressure and consequently 
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Dimensionless wavenumber q 

Figure 3. Solutions to the dispersion relation between the 
squared growth rate s'^ and dimensionless wavenumber q for a 
stellar disk. Solid curves correspond to Qg = 0.3 and the dashed 
curves are for Qg = 0.5. The dotted lines correspond to Newtonian 
gravity. 


the sound speed for this system is zero. However, because 
of the epicycle motion at any given patch in the disk there 
are stars from different parts of the perturbed disk. This 
cause an important cancellation effect (Toomre 1964). In 
fact, it turns out that the stellar disk can be described 
by the Jeans equations, which their mathematical form 
is rather similar to the hydrodynamic equations m and 
®. However, although equation (l20l) does not change, 
equation (EU is replaced with (see Roshan & Abbassi 
(2014b) for more detail) 

Vra = -^- k^aJ^ (37) 

where J- is the reduction factor given by (Binney & 
Tremaine 2008) 

-I _ 2 /“TT 

= sin sr sinrdr (38) 

sinTTs Jq 

where s = —is' = {uj — mVL)/n and q = kurj k, and ar 
are the dimensionless wavenumber and the radial veloc¬ 
ity dispersion, respectively. Since the reduction factor 
does not depend on the form of Poisson equation (Roshan 
& Abbassi 2014b), in MOG it possess the same general 
form and exhibit similar dependencies on dimensionless 
frequency and wave number as it has for Newtonian grav¬ 
ity. Now, substituting (1^ into equation (12(1 . we get the 
following dispersion relation 

(w - mVif = - 27rGEo|fc|[l + a - y=L=]T{s,q‘^) 

y/k^ + Mo 

(39) 

at the limit a —i 0 or /tq —>■ 0, the standard dispersion 
relation is recovered. Again, we restrict ourselves to ax- 
isymmetric perturbations. The system is stable against 
all axisymmetric modes if equation (1391) does not have a 


solution with negative In order to check this expec¬ 
tation, we rewrite dispersion relation (1391) as follows 


27 r|g| / _ a\q\ \ T{-is',q^) 

3.36Qs y \/q^ -\- j 1 + 


(40) 


in this case /3 = k. It is straightforward to ver¬ 

ify that the maximum value of the RHS of (14(1 occurs 
at s' = 0. Therefore, if the RHS of equation (HU for 
s' = 0 is smaller than 1, then there is no perturbation 
with negative which satisfies the dispersion relation 
( 001 ). Consequently, the disk will be stable against all 
perturbations. Therefore, the stability criterion reads 


Qs > 


27r|g| 

3.36 


I + a — 


a\q\ 




J-(0,q2) (41) 


This criterion can be compared with the stability crite¬ 
rion for a fluid disk, i.e. equation (031) ■ The criterion 
gB) can be replaced with a more simplified form 



(42) 


where /q is the modified Bessel function of order 0. In 
deriving (I42|) . we have used the following identity 

T{Q,q^) = ^{l-e-‘^"lo{q^)) (43) 

Equation (HU) is the main result of this section. For a 
given perturbation with wavenumber g, if we know the 
parameters a and /3 then we will easily calculate the re¬ 
quired value of Qa for stability of that perturbation. On 
the other hand, if we could find the maximum value of 
the RHS of equation (1421) , then we will obtain a criterion 
which guarantees the stability against all wavelengths at 
every location on the disk. In Table [TJ the maximum 
value of the RHS have been calculated for various values 
of a and /3. The tabulated data show that the required 
values of Qs for stabilizing the disk is larger than those 
for Newtonian case. As we mentioned before, this fact is 
also the case for a fluid disk. 

In order to compare the growth rates in fluid and stel¬ 
lar disks, we have plotted the stellar disk dispersion re¬ 
lation (HUl) in Figure [H This figure can be compared 
with Figure [H The solid curves correspond to Qs = 0.3 
and P = 0.04 to 0.2. The dashed curves correspond to 
Qs = 0.5 and P = 0.04 to 0.2. Furthermore, the dotted 
lines are the corresponding Newtonian dispersion rela¬ 
tions for the above mentioned values of Qs- As in the case 
of a fluid disk, the growth rates increase with decreasing 
parameter Qs- Also larger P leads to larger growth rate 
for large wavelengths. In fact, from Figure [3] it is evident 
that for small wavelengths there is no difference between 
the growth rates in MOG and Newtonian gravity. This 
point is also completely the same as for the fluid disk. 

It is also instructive to plot the boundary between the 
stable and unstable modes. As we mentioned before, this 
curves can be specified by the requirement uP = 0. In 
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Figure 4. The parameter /3 for stellar and gaseous components 
of some spiral galaxies. For gaseous component /3 increases with 
radius. For the stellar component ^ first increases and for large 
radius declines. 


Figure 5. Ratio of the critical surface density of Newtonian grav¬ 
ity and MOG. The red circles shows the ratio for the stellar com¬ 
ponent and the green circles correspond to the gaseous component. 


this case, the dispersion relation (1401) can be rewritten as 
f{Qs,y)y . , a 

i-ef(Q^’y)lo(f(Qs,y)) I ^A+^V7(Q^T) 

(44) 

where f{Qs,y) = • The curves corresponding 

to (H41) have been plotted in Figure [2] with solid lines for 
different values of /3 (0.02 to 0.2). It is obvious that sim¬ 
ilar to the fluid case, the range of instability is stretched 
to larger wavelength. The largest unstable wavelength is 
A = 10.23Acrit (note that the we have chosen a = 9.23). 
However as we discussed in the previous section, this ob¬ 
vious deviation from the standard case does not neces¬ 
sarily lead to a significant difference between MOG and 
Newtonian gravity. 


5. DISCUSSION 


just for illustration, and our results in this section are 
not restricted to these galaxies. For every galaxy in this 
catalog, the rotation curve v{r), epicycle frequency K{r) 
and the velocity dispersion ar(r) and Vs are known. So 
we can easily calculate /3 with respect to r for both the 
stellar and gaseous component of the galaxies. 

Before computing /3 from observational data, it would 
be useful to make an order of magnitude estimation. It 
is worth mentioning that, at large radii in spiral galax¬ 
ies the rotation curve is flat, i.e. v{r) ^ iiflat where 
Liflat is constant. Therefor, r2(r) ^ 2^, and conse¬ 
quently K(r) = For gaseous component of the 

galaxies, Vs is almost constant with respect to r and is 
Ts ~ 11 km s“^. In fact we have assumed that the sound 
speed is equal to the gas velocity dispersion Ugas, see 
Leroy et al. (2008). Thus at large radii /3 rises linearly 
with radius as 


/3(r) 


ypVs 

V^Uflat 


(45) 


As we showed in sections [3] and 01 the deviation be¬ 
tween the stability criteria of Newtonian gravity and 
MOG is highly controlled by the magnitude of param¬ 
eter /3. Therefore, in order to compare these theories, 
it is essential to find the magnitude of this parameter 
with respect to radius from observational data in some 
real spiral galaxies. To do so, we use a sample of six 
galaxies of THINGS catalog of spiral galaxies (Leroy et 
al. 2008). We have chosen this sub-sample of galaxies 


It means that the magnitude of /3 is grater near the edge 
of the disk galaxy and consequently higher values for 
Qg is required to overcome the gravitational instability 
there. An immediate and important consequence follows 
from this equation. Near the edge of the galactic disk the 
difference between Newtonian gravity and MOG might 
be significant. Therefor, the outer disks of spiral galax¬ 
ies are exposed to instability more than the inner disks. 
Gonsequently, there might be some regions near the edge 
























of the disk which can be stable in Newtonian gravity but 
unstable in the context of MOG. Such a deviation be¬ 
tween these theories at large radii is expected because 
the free parameter ^ 24 kpc is a large length scale 
compared to the characteristic size of the galaxies. A 
a consequence, taking into account the presence of the 
Yukawa term in this theory, we know that the effects of 
MOG appear at large radii. 

From TablelU it is evident that the stability criterion is 
significantly different from the standard case if /3 > 0.2. 
Using equation (145 1) the requirement /3 > 0.2 reads 

r > 6.7 kpc (46) 

If we set an average value ^ 10, then deviation 

between Newtonian dynamics and MOG will appear at 
r > 67 kpc. On the other hand, almost for all spiral 
galaxies, the characteristic optical radius is smaller than 
67 kpc. Therefore, this quick estimation shows that it is 
unlikely to find a galaxy where the predictions of MOG 
are substantially different from that of Newtonian grav- 
ity. 

However, let us do a more careful analysis by finding 
function j3(r) for our sample. The rotation curve v{r) for 
the spiral galaxies of THINGS catalog can be written as 
(Boissier et al. 2003) 

—r 
^flat 

where Wflat and Iflat are free parameters that refer to the 
velocity at which the rotation curve is flat and the length 
scale over which it approaches this rotation speed, re¬ 
spectively. These parameters are not universal and are 
different from galaxy to galaxy. For these parameters we 
use Table 4 of Leroy et al. (2008). Also, using the ro¬ 
tation speed v(r) the epicycle frequency can be written 
as 

K{r) = lAl^VTTb (48) 

r 

where b = dlnv{r)/dlnr. Furthermore, the radial veloc¬ 
ity dispersion is given by (van der Kruit 1988; Leroy et 
al. 2008 ) 

ar(r) = 1.55-\/ GLE*(r) (49) 

where S* is the stellar surface density and is the stellar 
scale length. We use Table 4 in Leroy et al. (2008) for 
For THINGS galaxies, E* has been measured with 
respect to radius in Leroy et al. (2008), for example see 
Table 7. Finally, using equations (H71) - (H^ . we have cal¬ 
culated fd for both the gaseous and stellar components 
of the sample. The result has been illustrated in Fig¬ 
ure |H The red circles show j3{r) for gaseous component 
and green circles correspond to the stellar component. 
As we have already expected, for the gaseous component 
P at large radii increases linearly with r. However, as 
illustrated in Figure IH although /3 increases with radius, 
it remains small even at the outer disks. This happens 
also for the stellar component. In other words, /3 remains 
small for stellar component. Of course, the behavior of 
P with respect to r for stellar component is completely 
different from that of the gaseous component. In fact. 



,{r) = 


I — exp 


MOG predicts that at outer disks the gaseous compo¬ 
nent must be more unstable against local perturbations 
rather than the stellar disk. However, as we mentioned 
above, this difference is not too big to be significant and 
detectable. 

In order to make the comparison between Newtonian 
gravity and MOG more clear, let us rewrite the stability 
criterion as follows 


E(r) < Ecrit(?') (50) 

In other words, a local axisymmetric perturbation at r 
is stable against gravitational collapse if the background 
matter density at that point is smaller than the critical 
density Ecrit(?')- Where, Ecrit(?') in Newtonian gravity 
for gaseous and stellar disk are respectively 

(Eerit)iV - , (Scrit)Ar - g (51) 

On the other hand using criteria m and dMl) , one can 
write the critical surface density in MOG as follows 

(Scrit)MOG = (^2) 

H{a,P) 

Where H > 1 and is a function of a and p. Note that 
for our purpose here, the exact form of the function H is 
not needed. Therefore, we can conclude that the critical 
surface density in the context of MOG is smaller than 
the Newtonian gravity. This makes a disk more unstable 
than it would be in the standard case. In order to see 
the difference between (Ecrit)MOG and (Ecrit)Af, we have 
plotted their ratio in Figure [5j Although for the gaseous 
component the ratio (Ecrit)Ar/(Ecrit)MOG increases lin¬ 
early, but it does not get significantly larger than 1. In 
fact this ratio can be approximated as 


(Ecrit)N 

(Scrit)MOG 


1 + /3 


(53) 


Since P remains small even at the outer disks, then one 
can conclude that the threshold matter density for grav¬ 
itational collapse in Newtonian gravity and MOG is al¬ 
most the same. Therefore, from this perspective, it seems 
like that there is no difference between these theories. It 
should be stressed that, we have done the same analysis 
for all galaxies in the THINGS catalog. And the result 
of this section is applicable to all of them. 


6. SUMMARY AND GONCLUSIONS 

In this paper, the local gravitational stability of galac¬ 
tic disks has been investigated in the context of MOG. 
The dispersion relation in the WKB approximation has 
been derived for both gaseous and stellar disks in the 
context of MOG. Using the relevant dispersion relations, 
we have found the local stability criteria for gaseous and 
stellar disks. 

Finally, we have used a sample of spiral galaxies of the 
THINGS catalog in order to compare these theories us¬ 
ing the current observational data. Our final results have 
been summarized in the Figures |4] and [5l In Figure [3 
we have compared the critical surface density over which 
the local perturbations collapse, in Newtonian gravity 
and MOG. For all galaxies of the sample, although the 
ratio of (Ecrit)Ar/(Elcrit)MOG increases with radius, but 
it does not differ substantially form unity. Therefore, at 
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least for all the spirals in the THINGS catalog, there is 
no significant difference between Newtonian gravity and 
MOG. It is also important noting that it could be more 
precise to find the stability criterion for a two-fluid disk 
(gas-fstars). However, taking into account the main re¬ 
sults of this paper, it is unlikely that a two-fluid analysis 
changes the results. 

Of course, our results in this paper do not mean that 
the dynamics of spiral galaxies are the same in Newto¬ 
nian gravity and MOG. Only the local perturbations fate 
is the same in these theories. It is worth remembering 
that MOG can explain the rotation curves of the spiral 
galaxies without using the dark matte halos. However, 
this can not be achieved in Newtonian gravity without in¬ 
voking dark matter particles. Also, one should note that 
the global stability analysis of the disk galaxies in MOG, 
in principle, can be different from the standard case, for 
example see Brandao et al. (2010) for a numerical study 
of disk galaxies in MOG. 

As a final remark, let us mention that there is a well- 
known instability in the dynamics of spiral galaxies called 
bar instability. In fact, in the early 1970s, the N-body 
simulations of self-gravitating disks showed that a cold 
rotationally supported disk is globally unstable to the 
formation of a non-axisymmetric system. In other words, 
a rotation dominated disk of particles will undergo a 
rapid transition to a bar-like pressure dominated system 
(for example see Hohl 1971). This result was puzzling 
because real galaxies like the Milky way does not show 
such a behavior. However, it is known that the presence 
of a dark matter halo can help to stabilize the disk, see 
Ostriker & Peebles (1973). Furthermore, the presence 
of a dense bulge-like mass component near the center of 
the disk can also stabilize a cold self-gravitating disk, see 
Sellwood (2014) for a comprehensive review of the sub¬ 
ject. Another approach to this problem can be provided 
via modified gravity theories. In alternative theories of 
dark matter there is no dark matter halo and the bar in¬ 
stability problem should be addressed without invoking 
dark matter. For example see Brada et al. (1999) and 
Tiret et al. (2007). In these papers the global stability of 
the disk galaxies has been investigated in MOND. As we 
mentioned before, MOG is an alternative theory of grav¬ 
ity presented for addressing the dark matter problem. 
Therefore, it would be interesting to study the bar insta¬ 
bility problem in this theory. The result of the current 
paper can be used in numerical study of disk galaxies in 
MOG. More specifically, in order to set the initial condi¬ 
tions of the particles in a N-body simulation for studying 
the global stability of the disk, the local stability crite¬ 
rion on the disk is required to be satisfied. 


We would like to thank the referee for useful comments. 
This work is supported by Ferdowsi University of Mash¬ 
had under Grant No. 100835 (25/05/1393). 
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